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We analyse an ambiguity in previous works on entanglement of fermionic fìelds in non-inertial 
frames. This ambiguity, related to the anticommutation properties of field operators, leads to 
nonunique results when computing entanglement measures for the same state. We show that the 
ambiguity disappears when we introduce detectors, which are in any case necessary as a means to 
probe the field entanglement. 
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I. INTRODUCTION 

We will discuss that for bipartite systcms of fermionic 
fìelds, anticommutation of fermionic operators naturally 
induces an ambiguity when defining the individuai basis 
of the observers. This ambiguity is related to the ordering 
criterion of the creation and annihilation operators and 
has gone unnoticed in ali previous literature on fermionic 
field entanglement in non-inertial frames [l|-[loj. 

Such ordering is only a mathematical convention since, 
obviously, physical states are the same regardless of how 
wc rearrange the fermionic operators. Therefore, one 
would expect that choosing a different convention should 
not changc the bchaviour of entanglement measures. We 
will show that using the formalism of ali the previous 
literature, the results do depend, in general, on this con- 
vention. 

However, to eventually measure field entanglement, we 
need to consider detectors coupled to the field. We argue 
below that the ambiguity is removed when we incorporate 
detectors to the field entanglement settings. 

Incidentally we show that the convention that has a 
physical meaning in terms of detectors entanglement is 
different to that used in the previous literature, where 
field entanglement (without detectors) was analysed in 
the Dirac case and Grassman fields beyond the single 
mode approximation P, . 

The paper is organised as follows: In scction we 
present and thoroughly analyse the ambiguity that ap- 
pears when defining bases from different ordcrings of fer- 
mionic operators. In section Hill we study the implica- 
tions of this ambiguity in the specific context of relativ- 
istic quantum Information, finding that previous liter- 
ature results are not independent of different sign con- 
ventions. In section HVl we show how for any physically 
meaningful experiment such ambiguity disappears. We 
suggest a physical criterion to choose the signature of 
the fermionic basis, based on meaningful outcomes of de- 
tectors experiments. Finally, conclusions are presented 
in section Ivi 



II. A TOY MODEL 

In this section we study the entanglement properties of 
a simple fermionic system to illustrate the peculiarities 
wc intend to show in more general settings. Namely, 



we consider a two-mode fermionic system with fermionic 
creation operators and 6^ , and a vacuum state 1 0) . The 
relevant Hilbert space T-L is thus four dimensionai. Any 
fermionic system in which there are only two possible 
one-particle states will fit to this quite general model. 
We denote the Hilbert space basis of our toy model as 



|00) = |0), |10) =at|0), |01) =òt|0) 
|11) =at5t|o). 



(1) 



By choosing this basis, we implicitly cndow the Hilbert 
space with a tensor product structure which allows us to 
regard it as two qubits, with the first label corresponding 
to one qubit and the second to the other. We can now 
perform quantum information and study entanglement 
properties in this system. 

Of course, if we make a nonlocal change of basis, the 
entanglement of the state may change [Ì4I. Due to the 
fermionic nature of our system, there is one such change 
of basis which arises naturally: we can interchange the 
positions of and in ([T]), or equivalently, we may 
change the last term by a sign. The new basis, which we 
denote by a prime, is 



|00)' = |00), |01)'=|01) 

|10)' = |10), |ll)' = òtat|0) = -|ll) 



(2) 



Thus, in this specific case, reversing the order of and 
leads naturally to a new basis, which only differs by 
a sign in one of its elements from ([T]). We usually think 
of the modcs and 6^ as cach one acting on one qubit. 
A naive cxpectation would be that the ordering betwcen 
them should be purcly convcntional and should not affect 
entanglement. Nevcrthclcss, a great deal changes from 
(HI) to ((21): The new basis endows the Hilbert space with 
a different tensor product structure. By this we mean the 
(quite obvious) fact that the set of elements of the Hibert 
space which are expressible in the form vi® V2^ with vi , 
1)2 of the form ci |0) -|- C2 |1), is different if we takc the 
basis (H]) or ([21). In the language of quantum information, 
the transformation from ([1]) to ([J) is not locai, and there- 
fore it cannot be expressed as the product of two unitary 
transformations, with each acting on one qubit. This 
means that, for any given state, entanglement properties 
will vary wildly in general, depending on which operator 
ordering is taken. 

We remark that this is true for any nonlocal change of 
basis in any bipartite quantum system. Indeed, we could 
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have imposed the change of basis from ([T]) to ([2]) in a 
bosonic system as well. Nevertheless, whereas in a bo- 
sonic setting such a change would have been arbitrary, in 
the fermionic case, due to the sign acquired by operator 
interchange, there is no a priori way to choose between 
these two bases to define our bipartite system, unless a 
specific set of obscrvables with a spccific operator ordcr- 
ing structurc is choscn (see section UV]) . 

As a concrete example, considcr the scparable state in 
basis di]). 



Ivi/) = 1(100) + 101) 



|io) + |ii)) 



V2 



(|0) + |1)) 



V2 



(|0) + |1)) 



(3) 



When expressed in basis ©, it reads 



l*) = 2(|00)' + |01)'+|10)'-|11)'), (4) 

which is a maximally entanglcd state as measurcd by 
the Von Neumann entropy of its reduced density matrix, 
which is 5 = 1. In fact, if we perform the change of basis 
IO') = ^ (IO) - |1)), |1') = ^ (IO) + |1)) in the second 
qubit space, the state |vl>) becomes the Beh state |'I'^). 

This shows that when studying quantum information 
in fermionic systcms it is fundamental to choose a spc- 
cific operator ordering. This phcnomcnon happcns in any 
fermionic system rcgardless of its origin, and in particu- 
lar in quantum ficld thcory of the Grassman scalar and 
Dirac fields, which is a point that has been neglected in 
ali previous works known to us [lUloj . 

So far, we have only described the naturai con- 
scqucnces of making a nonlocal change of basis in a bi- 
partite system. It could be argued that as long as we stick 
to onc basis, albeit arbitrary, we will obtain well-defined 
rcsults. However, as we show below, entanglement can 
also change when partial traces are taken, even if we do 
not change the operator ordering in the remaining state. 
Thus, the 'ordering convention' in the unobserved degrees 
of freedom may scem to affect the entanglement behavior 
of the rest of the system. 

Consider a tripartite system. If for any reason we do 
not have access to the third subsystem, i.e. ali of the 
observables that we consider belong only to the first and 
second subsystems, then we have to sum over ali of the 
States of the unobserved Hilbert space. Formally, this is 
done by taking a partial trace over the third subsystem. 
After tracing, we end up with a bipartite state which is 
gcncrally nonpure. 

In the absence of a physical critcrion to adhcrc to a 
particular 'operator ordering', there is no reason to re- 
strict ourselves to a fixed Hilbert space basis. Therefore, 
we will choose the naturai basis associated with each op- 
erator ordering in the same way as we did in equations 
di]) and dU for the orderings a^b^ and b^a\ respectively . 
From now on, when we speak about 'different orderings', 
we are implying that our sign convention when defining 
the Fock spacc basis is adapted to that ordering. 



We will now show that even if the untraced modes 
do not change their relative ordering, entanglement may 
change depending on the 'ordering convention' in the 
third mode. In other words, entanglement is depend- 
ent on the relative position of traced out modes, which 
are no longer present after the partial trace is taken. 

We will show this phcnomcnon by minimally modifying 
our toy model: Let us add a third mode created by c^, 
and considcr the state 

l'J'> = ^(|100) + 1010) -1-1101) + 1011)) (5) 

where the basis is defined using the ordering (a^ò^c^), 
so as to have |111) = a^b^c^ |0). Suppose that the c' 
modes are not observed and thus results are obtained 
after tracing over them. This means tracing over the 
third label in dS])- To study entanglement between the 
first two subsystcms wc havc to take a partial trace over 
the third onc. Wc then study the negativity [l^ of the 
reduced mixed state, which would contain only and b^ 
excitations. The negativity for state dS]) after tracing out 
the c' mode is i, so the state is entanglcd. 

Of course, we expect that interchanging the ordering of 
and will result in different negativities. This is so, 
since entanglement is dependent on operator ordering as 
we showed above. However, it is also true that even if we 
do not interchange these operators, and only permute 
with any one of them, entanglement on the reduced state 
changes. This is remarkable because disappears after 
taking the partial trace. As a matter of fact, if we take 
the ordering (a^c'b^), so as to have |111) = a^c^ò^ |0), the 
last term in dS]) changes sign and the reduced state is 



P = 2(|10)(10| + |01)(01|) 



(6) 



which is a complctely uncntanglcd state. The negativity 
is thus zero. 

A situation of the kind described above appears in the 
context of quantum ficld thcory in non-inertial frames 
and curved spacctimes; see section HIT] and [il, I3|- Spe- 
cifically, in certain cases (uniformly accelerated observ- 
ers, entanglement in the presence of black holes, etc), 
there are regions of spacetime which are causally discon- 
nected from the observer's world-line. Therefore, modes 
living inside them are not allowed to communicate with 
the obscrver. One is therefore forced to trace them out 
[il, [l3|. In the context of fermionic fields, the relative po- 
sitions of these modes will affect entanglement, as above. 
Thus the entanglement in the ficld state dcpcnds on the 
sign criterium over modes which do not have a causai 
connection with the observer. 

It is obvious that physical observations cannot be 
changed by an operator ordering criterium. Indeed, a 
specific ordering is imposed by the physical nature of 
experiments as we show in section IIVI However, this 
fact has been ovcrlooked in previous studies in relativ- 
istic quantum information, assuming somchow that an 
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arbitrary convention can be taken to derive general res- 
ulta. Although this is true in some cases, we will show 
that remarkable differences appear when we introduce a 
physical criterion to select the basis. 

We conclude this section with the rcmark that there 
does not seem to be any way to construct an ordering- 
indcpendent entanglement measure for fermionic states. 
Any reasonable entanglement measure should vanish for 
separable states and achieve its maximum at maximally 
entangled states. But as our examples show, separable 
states in one basis need not be separable in another, and 
can in fact be maximally entangled. It is therefore im- 
possible to define a reasonable entanglement measure for 
these states. 



III. OPERATOR ORDERING IN 
RELATIVISTIC QUANTUM FIELD THEORY 

The effects of operator ordering will now be discussed 
for the kind of states that appear in relativistic quantum 
Information. Specifically, we will consider a 1+1 fiat 
spacetime in which there is an observer, Rob, who undcr- 
goes uniform proper acceleration, and another observer, 
Alice, who stays inertial. Both look at the same state of 
a fermionic field of spin s, but while Alice uses a basis 
of Minkowski modes (piane waves in the massless case) 
for the description of her part of the field state, Rob uses 
the so-called Rindlcr modes Q. These modes are the 
naturai candidates to describe the state of a quantum 
field from the viewpoint of an accelerated observer, since 
a uniformly accelerated detector would couple to them. 
Due to the nature of the change of basis between Rind- 
ler and Minkowski modes (give n by the so-called Bogoli- 
ubov cocfhcients; see 

dMlìll), which mixes Minkowski 
creation and annihilation operators, the Minkowski and 
Rindler vacuums are not the same. This is the origin of 
the celebrated Unruh effect . 

Rindler modes are named after Rindler coordinates 
(?7,r), which describe a family of uniformly accelerated 
observers in two causally disconnccted patches of space- 
time, labelled as regions I and II in Figure [T] The ob- 
servers are right-moving in region I and left-moving in 
region II. An observer foUowing an orbit of Constant rj 
is uniformly accelerated and sees the boundary of region 
I as an evcnt horizon. Likewise, the same happens for 
an observer foUowing an orbit of Constant rj in region II. 
Rindler modes are monochromatic solutions of the field 
equation in Rindler coordinates, therefore they only have 
support inside either region I or II. This means that for 
each particle species, frequency w and spin-z component 
we have two annihilation operators, one for each region. 
The operator c^,s.i will corrcspond to the Rindler mode 
of frequency lo and spin a in region I, and d^^s,i will cor- 
respond to its antiparticlc, with the same considerations 
applying to region IL 

The existence of these two kinds of modes implies 
that, from the accelerated observer viewpoint, the Hil- 



Rob 



/■■ I 



Alice 



Figure 1. Minkowski spacetime diagram showing the world 
lines of an inertial observer Alice, and one uniformly accel- 
erated observer moving hyperbolically in region I. Note that 
regions I and II are causally disconnected from each other. 



bert space factorises as Hi (E) Hu- Any physical acceler- 
ated observer such as Rob must remain in either region I 
or II, since both are causally disconnected. In quantum- 
mechanical terms, this means that we must trace out the 
part of the state outside the observer's region, since it 
will not be observed. Customarily we take the observer 
Rob to be in region I, and so henceforth we will trace out 
-Hii. 

There is another set of modes which can be constructed 
by taking simple linear combinations of Rindler modes 



^l,a.R = cosr^^cl'^^^ j - sinr^(i„,_cr,ii, 



Cia.L cosr^c;';, jj - sinr„d^,_^,i, (7) 



where tanrj^ = ^-TTuic/a ^^^-^^ ^ indicates Rindler fre- 
quency. These are the so-called Unruh modes. The 
subscripts L and R stand for 'left' and 'right' modes, 
which are related to each other by a reversai of regions I 
and II. These modes have the particularity of being lin- 
ear combinations of purely positive-frequency Minkowski 
modes. In terms of creation and annihilation operators, 
this means that ([7]) can be rewritten as a linear combin- 
ation of only Minkowski creation operators. This trivi- 
ally implies that the Minkowski and Unruh vacua are the 
same. The Minkowski vacuum can then be written as a 
product of the vacua for each Unruh mode. 



IO) 



M 



IO) 



(8) 
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where each |0)u is annihilated by C^j.^-.r and C^^^o-.l for 
ali a. 

In this work (as in prcvious litcrature [H, [Ì3l): the ac- 
cclcratcd observer Rob watches a single Unruh mode of 
the field. Given (|5]), we only need to perform a detailed 
study of one of the factors in the tensor product, and 
hence from now on we shall drop the frequency index 
in ali operators. We will consider arbitrary Unruh ex- 
citations instead of restricting ourselves to the so-callcd 
single mode approximation [9| 



|lU = <s,ulO>u 



R + (ii^cl ^ A |o)u , 



(9) 



The case qn = 1 corresponds to the choice previously 
known as the single mode approximation [isl [l9j . 

Now, the Unruh vacuum and excitations are expressed 
in a simple way in the Rindler basis, due to the trivial 
form of ([7]). Note that no considerations whatsoever of 
operator ordering have been assumed so far: To fix a con- 
vention, we will assume Alice's mode appears leftmost in 
ali expressions. It is also assumed that the mode watched 
by the inertial observer Alice possesses a negligible over- 
lap with the Unruh mode, so that we can consider both 
subsystcms as noninteracting. 

For concreteness, we will be interested in states of the 
form 

Ivi/) = P |0)a (4 |0)u) + Q |1)a (4 |0)u), (10) 



= 1, 



where and B\j are arbitrary linear combinations of 
products of Unruh modes C„^\j so that Rob's state is in 
general not a qubit. Note that the second part of the 
state will have to be expressed in the Rindler basis. 

In order to express |0)y in the Rindler basis, we will 
factor the vacuum as |0)u = |0)p (8) |0) j with Ccr,R |0)j^ = 
and Ccr,L |0)l = for ali a as in Q. This factorisa- 
tion alrcady entails an assumption on the ordering of 
the Rindler operators. Namely, the operators ^ j and 
should always appear to the left of the operators 



" LJ, — (7,11 



c, , ^ TT and d,, „ t. 

We will now choose a specific ordering for both the 
left and right sectors as follows: within each sector, ali 
region I operators shall appear before region II operat- 
ors and. within each region, operators will be ordered by 
their spin-z component, with the highest vaine appearing 
leftmost. We can associate a binary number to each ele- 
ment of the basis. In this fashion, the vacuum state for 
each sector would be indexed by a sequence of 2 • (2s -|- 1) 
zeros. The one-particle excitations would have a 1 in the 
fcth position, and so on. For example, for s = ^, we 
would have 



uni) 



^ ^t-l'^lj'^TJl'^LlI l'^)Rindlcr' 



|llll)L = 4.i4.i4,n4iilO) 



Rindler * 



(11) 



This notation coincides with that used for the Grass- 
man field in prcvious works Q and has the advantage 
that it is easy to find and write the vacuum state for 
arbitrary spin, as we will show below. 

In the litcrature for spin 1/2 fields, [H the common 
notation groups the spin-up and spin-down operators for 
each region and particle species. To translate from the 
generic notations introduced above to the short one for 
Dirac fields, we take the indices in groups of two and 
make the replacements 00 — >■ 0, 01 — >-4,, 10 and 
11 p. The p stands for 'pair'. For instance, we have 

iioii)j,|oiio)L = itp;t) 

We now can calculate the right vacuum by assuming 
the foUowing ansatz, in which again we have introduced 
new notation: 



|0)i 



(12) 



Here, a is a binary number with length 2s -|- 1, and 'R-ipì) 
is the binary number obtained by reversing the ordering 
of the figures of a [e.g. if a = 01 TZia) = 10]. The sum 
is extended to ali 2^"+-'^ such chains. The chain olR,(oÌ) is 
obtained by concatenating them so that a is associated 
with modes in region I and 72.(a) with modes in region 
IL The ansatz (IT^ is of the form of a squeezed vacuum 
state. 

We shall define a few operations and notation on these 
chains: x(q^i fc) will denote the number of l's in a before 
its fctli position. x(q^) will simply be the sum of ali the 
digits in a. If a is a chain with a in its (2s -I- 2 — fc)th 
position, then a -|- 2*^ will be the chain with a 1 in its fcth 
position and the same digits as a elsewhere, mimicking 
the addition of binary numbers. Finally, Su will be the 
set of ali a with a in its (2s -I- 2 — A:)th position. 

In ordcr to determine the coefficients Xa wc impose 
C'o-.R |0)p{ = for ali a and therefore get 







OL 

E 
E 



X^I^COS 7'^CcrJ |a Ti-icì)) — 



^-cr,II 



|a ??.(«)) j 



Xa cosr^jCcr,! la ^.{a))^ 
Xa sinr^d^^ jj \a TZ{a))^ 



(-l)x(»^-)xa+2^ cosr^ \a 7^(a + 2"))j 



= E 

(„l)x(a)+x(K(a),2s+2 



' Xa sm r. 



aes. 



(-1)> 



+2" 



a TZ{a -f 
Xa sin r^) 



2^))i 



|« 7^(a + 2-))j, 



(13) 



In passing to the last line, we used the property 
X(7^(a),2s + 2 - (j) = x{a) - x(a, cr) for ali a e Sa- 
By equating each coefhcient to zero wc obtain the recur- 
rence relations 



— tanr^ • Xq,, 



(14) 
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which can be solved iteratively by taking, e.g, xo,,,o = 1. 
Then, Xq. = tan , and upon imposing normalisation, 
we obtain for the right vacuum 

|0)r = ^(cosr^)2^+i-'^(")(sinr^)'^(") \a n{a))^ (15) 

a 

The left vacuum is straightforward to obtain from (|15p 
by noting that the ordering in the left sector intcrchanges 
particles and antiparticles. The left vacuum thcrefore is 

|0)L=^(-l)'^("Hcosr„)'"+i"^("Hsinr„)>^(") \an{a))^ 

a 

(16) 

We remark again that in deriving ([TS|) and ((T5)) we relied 
on a very specific operator ordering. As seen in sectionllll 
entanglement measures will not remain invariant in gen- 
eral when this ordering is changed. This is a point so far 
neglected in the literature.We will now study the changes 
in negativity for the state ([TU| and different values of s. 

A. Grassmann fìeld 

We shall now consider the bchaviour of entanglement 
under different operator orderings in the Grassman scalar 
field, which is an anticommuting scalar field with only 
one degree of freedom. This field has been extremely 
useful to study the gen eral features of entanglement in 
fcrmionic fields [Ì1,|1IÌ3- 

Under these circumstances, we can use ([TS]) and (fT6|) 
together with \0)^ ~ |0)r ® |0)l obtain the Unruh va- 
cuum in tcrms of the Rindler vacuum, using the operator 
ordering cj djj cjj: 

|0)u = cos^ 10000) — sinr;^ cosT;^ |0011) 

-hsinr^cosr,^ |1100) - sinV^ |1111) . (17) 

The one-particle excitations are obtained as 

|1)u = (9R^1 + 9lCÌ)|0u) 

= QR [cosr^ 11000) - sinr„ |1011)] 

+ QL [sinr^ |1101) + cosr„ |0001)] (18) 

and the general state pU)) takes the simple form 

1*) = P |0)a [ai |0)u + h |l)u] + Q |1)a [«2 |0)u 

+62|l)u], |a,|2-f|ò,|2^1 fori = 1,2. (19) 

We now compute the density matrix for the state, and 
then trace out region II operators for the reasons de- 
scribed above. Finally we compute the negativity as an 
entanglement measure for the Alice-Rob system. For a 
general state and values of qr, entanglement can be cre- 
ated as shown in [20j , in contrast with the common con- 
ception that the Unruh effect should effcctively act as a 
thermal bath. 



However, we could do this in a different operator or- 
dering. For instance, we could rearrange the region II 
operators in any way we want while leaving the relative 
positions of the Alice mode operator and Rob's region I 
operators unchanged. A naive expectation would be that 
since region II modcs live in a different causai patch of 
spacetime, thcir position should be irrelevant when con- 
sidering entanglement. 

Nevertheless, this is not the case. The expression for 
the state (|19p in any operator ordering is readily obtained 
from ((T7)) . (US]), and ((T^ by permuting the operators in 
each term of the superposition, and taking into account 
the corresponding signs. Then, when region II modes are 
traced, the phenomenon discussed at the end of sectionllll 
appears: There are indeed different negativity behaviours 
corresponding to different operator orderings, as shown in 
Figl2] Note that in the Grassman case, only two different 
behaviours appear for a general state. This can be traced 
back to the fact that by performing some change of basis 
we can get rid of ali the signs in the partial transpose of 
the reduced density matrix, except for one. 



0.26 




0.16- 
0.14- 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 
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0) 

Figure 2. (Colour online) Negativity as a function of rui for 
QR = I/V2, the state ([T9| with P = 1/^2, «2 = òi = 0, 
ai = 62 = 1 and operator orderings cj dfj c?J cjj (red dashed 
curve) and djj cjj (blue solid curve). Note that in the first 
case, entanglement is a monotonically decreasing function of 
Tu and thus acceleration, whereas in the second, entanglement 
is generated for r„ > ^. 

There is always some negativity that survives in the 
infinite acceleration limit, in concordance with previous 
results [H-i, i, [H HH, This is true for ali the states 
and fcrmionic fields we have considered. The survival of 
entanglement at infinite acceleration is therefore a phys- 
ical phenomenon and not an artifact of choosing a spe- 
cific ordering. Moreover, there is evidence that the phe- 
nomenon is related to a tradeoff of entanglement between 
particle and antiparticle sectors [lTI |. 

Remarkably, the state shows entanglement regencra- 
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tion in one ordering which is absent in the other. Our 
results in section |IT] imply that any other entanglement 
measure will present this ambiguity. What is needed is 
a physical criterion for choosing a specific ordering; see 
section IIVI 

Finally, we note that if we restrict ourselves to = 1 , 
the case prcviously known in the htcrature as the single 
mode approximation, there is only one negativity beha- 
viour. This is a quirk of the Grassman field and it is not 
true for general spin, as we shall see below. 

B. Dirac fìeld 

The Dirac field presents a much richer zoo of negativity 
behaviours, as could be expectcd from the doubling in the 
number of Rindler operators from the Grassman field. 
The general state (|T0|) now takes the form 

1*) = P |0)a [ai |0)u + fei |t)u + Ci |;)u + di b) J 

+ Q |1)a [a2 |0)u + b2 |t)u + C2 |;)u + d2 |p)u] , 

(20) 

where 

|a,p + |ò,|2 + |c,|2 + = 1 for i = 1,2, (21) 

and we bave defined \p)^ = Cj ^C] ^ |0)u. The vacuum 
and excitations can be computed from and ((T5)) in 
the same way as in the Grassman scalar case. The im- 

plicit ordering is cj ^ cj j dj ^ djji 4.1 4,i 4." 4,ii- 
order to make the results more readable, we shall use the 
abbreviated notation for the Dirac field defined above. 
The translation from the notation in which (fTS)) and p6)) 
are written to the new one is straightforward. With these 
convcntions, the vacuum, excitation, and pair terms are 

1 0) u = cos-^ 1 0000) - cos^ sin ( 1 00 ti) + 1 00 ;t) ) 
+ cos^ sin^ \00pp) 
+ cos3r„sinr„(|n00) + |;tOO)) 

- sin^ cos2 ruinu) + itut) + ìim) + imt)) 

+ cos sin^^ r^{\tlpp) + \l^pp))+ sin^ \pppp) 

- cosr„ sin^ r„(|pptl) + |pplt)) 

+ cos^ sin^ |ppOO) , (22) 

k)u = 9l[cos^ |000cr)+cos2 sinr<^(|t|Ocr) + lltOcr)) 
+ cos sin^ r^^ bpOa) ) 
+ sgn((T) • (cos^ sinr^ |GOctp) 
+ cosr^ siii^ r^(ltlo-p) + ap)) + sin^ r„ \ppop))] 
+ gR[cos3r^ IcrOOO) — cos rtjSÌnr(^(|crOti) + |(TO|t)) 
+ cos sin^ I^Opp) 
+ sgn(o-) • (cos^ sinr^^ IpcrOO) 

+ sin^ r„ \papp) - cosr^ sin^ r^dp^ti) + \p(^ V\)))], 

(23) 



b)u = 9r (cos^ bOOO) - sinr„ cosr„(|pOti) + IpOit)) 
+ sin^r^ IpOpp)) 

+ ql (cos^ r„ |000p) + sinr„ cosr^dt^Op) + \i^Qp)) 
+ SÌV? \ppOp)) 

+ gR9L (cos^ ItOOl) - cosr^sinr^ |tO|p) 

+ sinr„ cosr^ btO|) - sin^ \p^ip) 

- cos^r^ IlOOt) - cosr^sinr^ IlOtp) 

+ sinr^ cosr^ b|Ot) + sin^ \p^p)) ■ (24) 

Using (|22p . (1231), and ([M)) . we can study negativity for 
an arbitrary state of the form (|20p and operator order- 
ing, by making adequate permutations. It is not straight- 
forward to compute explicitly how many negativity be- 
haviours can arisc, sincc there are 8! = 40320 possible 
ordcrings. Many of these are cquivalent rcgarding entan- 
glement: For instance, any permutation which involves 
only transpositions of operators of the same region con- 
stitutes a locai unitary and thus leaves entanglement un- 
changed. For some special states, there can even be more 
symmetry: If (?r = 1, then ali the 'left' operators are 
irrelevant, meaning that their position in the ordering 
does not affect negativity. The converse is true when 
= 1- If the state only carries spin-up excitations, ali 
the spin-down particle operators and spiii-up antiparticle 
operators are irrelevant. 

Nevertheless, a survey of ali the 8! orderings was 
carried out numerically. The results vary wildly de- 
pending on the degree of symmetry of the state. For 
the singlet state, (1^01) with ai = ci = di = 02 = 
62 = ^2 = 0, òi = C2 = 1 and P ~ 1/a/2 there 
are only six negativity behaviours, as depicted in fig |3l 
At first negativity rises steadily, then reaches a max- 
imum after which it decrcascs monotonically. For the 
canonical ordering implicit throughout this work, neg- 
ativity tends to its inertial value of ^ from above as 
r,^ — >■ 00, never decreasing past this value. For the or- 
dering cJ j d| jj 4,1 4,11 4,1 4;!! 4,1 4,11 negativity de- 
creases far more quickly, reaching only 0.1398 in the in- 
finite acceleration limit. AH the other behaviours are 
intermediate between these two, and are qualitatively 
similar to the second. The six classes are not evenly 
populated: The uppermost one contains 4032 different 
operator orderings and is the least abundant, while the 
bottommost one contains 9408 operator orderings and is 
the most abundant. 

If wc consider a general state excluding pair terms, that 
is, a state of the form ((20|) with di = o?2 = 0, the num- 
ber of negativity classes rises to 64, which suggests that, 
in analogy with the Grassman case, the choice between 
one behaviour and anothcr is govcrncd by six signs. This 
is, however, disproved by the differing number of order- 
ings which show each behaviour: If the choice depended 
only on six signs, then ali classes should be equally pop- 
ulated. We also note that for (/r = 1, the 64 behaviours 
coUapse into only two, which is not surprising given the 
high degree of symmetry present in this case. 
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Figure 3. Negativity as a function of r^^ for the Dirac singlet 
state when = 1/ -\/2 and different operator orderings. The 
general behaviour is similar in ali such cases: There is a slight 
entanglement creation wliich is then destroyed. For every 
ordering, some entanglement survives ai Vui = ^ . 



infrequent behaviours. There are also families of differ- 
ent negativity behaviours with the same population. The 
naturai explanation for these phenomena is that each of 
these families correspond to a set of orderings with the 
same symmetries. 



X 10 




Different behaviour classes 



Figure 4. Population for some negativity classes in the s = | 
case, including the most abundant ones. Note that the points 
group in 'steps' of classes with the same negativity, and the 
two dominant classes are well above the others. 



Finally, considering ()20|) with arbitrary coefficients res- 
ults in 778 classes, which lack any exploitable structure. 
There are two effects determining the number and pop- 
ulation of the classes: On one side, we bave the different 
signs which can arise for each term of the superposition 
in (PU|) for each operator ordering. On the other, the 
density matrix (and thus the partial transpose used in 
computing negativity) is zero in most of its entries. This 
makes some sign changes trivial, which could otherwise 
bave resulted in a different negativity. 

C. Higher spin 

A similar treatment to those of the Grassman and 
Dirac fields can be performed for fermionic ficlds of half- 
integer spin. The complcxity of the task grows quickly: 
There are [4- (2s-|-l)]! possible orderings for a field of spin 
s. This means that already for s = |, for which there are 
2.1 • 10^'^ orderings, the use of Monte Carlo algorithms is 
essential. 

We considered the particular instance of pU)) given by 
P = I/V2 and 

^" = ^+iu + ^+i,u' ^u = ^-i,u + ^-iu (25) 

This state is equivalcnt to a Dirac singlet state by making 
a simple change of spin basis, but the number of negativ- 
ity behaviours is much higher because of the greater num- 
ber of operators (and thus possible permutations). We 
found ss 1.4 • 10^ different behaviours for negativity. As 
can be seen in Fig|31 there are two dominant behaviours 
much more frequcnt than the others, and a plethora of 



IV. CHOOSING A PHYSICAL ORDERING 

We bave seen in scction lTTll that there can be an enorm- 
ous number of different negativity behaviours for each 
different operator ordering that shuffles region II operat- 
ors. This leaves us in an uncomfortable position, since 
operator orderings are purely conventional and therefore 
a priori equivalent. But if ali the orderings seem to be on 
an equal footing, it is only because we are lacking a phys- 
ical criterion for which class of them should be chosen. 
We will find such a criterion in this section. 

First, note that ali observables such as expected values 
or transition probabilities are obviously manifestly inde- 
pendcnt of operator ordering: One can always use the an- 
ticommutation relations to transform any matrix element 
into a superposition of terms proportional to (0 |0) = 1, 
and this transformation is unique. Although negativity 
is not an observable, it can be experimentally determ- 
incd (for instance by means of state tomography) , which 
means that it ought to be expressible in terms of expect- 
ation values. Nevertheless, it is not the negativity of the 
field that which can be measured, but that of the de- 
tector. Indeed, for any particular detector-field model, 
we look at the state of the coupled system, usually long 
after the interaction has been turned off, and then trace 
out the 'field' subsystem. The negativity of the detector's 
entangled state can then be experimentally determined. 
A field is never really observed; it interacts with a de- 
tector which is the only experimentally accessible system. 

This means that in order to determine a suitable order- 
ing we must study the form of the interaction bctween 
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the detector and the field. A physical detector for fer- 
mionic fields has not yet been proposed. Nevertheless, 
there is an obvious causahty condition that any physical 
detector shoiild fulfiU: The interaction Hamihonian for a 
detector moving in region I of Rindler spacetime cannot 
contain any region II operators, and vice versa. 

We seek to find the conditions in which ah of the ficld 
entanglement is physicaUy meaningful, i.e., it can be ob- 
served. Without a fermionic detector model at our dis- 
posai, we cannot address this question directly; it may 
well be the case that no physical detector can ever ac- 
quire ali of the field correlations studied bere and else- 
where. In this case, we would be justified in labeling this 
field entanglement as unphysical. 

However, it may be the case that some particular de- 
tector can be imprinted with ali of the field entangle- 
ment. If so, there will be a particular operator ordering, 
dependent on the exact nature of the field-detector in- 
teraction, in which no anticommutation signs will arise 
in any computation (simply take the operator ordering 
to be the reversc of the fermionic operator ordering in 
the interaction Hamiltonian) . Since the interaction only 
has region I fermionic operators, any such ordering must 
bave ali region II operators rightmost. Permutation of 
the operators of the same region constitutes a locai unit- 
ary transformation, which docs not changc entanglement, 
and therefore any two such ordcrings will show the same 
entanglement. 

We therefore bave found a necessary condition for field 
entanglement to be physical. This condition selects one 
and only one entanglement behaviour class and therefore 
settles the question of what ordering should be chosen to 
computo fermionic entanglement. We also point out that 
the rcsulting negativity in states such as (fT9)) is diffcrent 
from the one reported in Q, which corresponds to a dif- 
ferent class of operator orderings. Namely, these works 
found the monotonically decreasing red dashed curve in 
Fig. [21 whereas the physicaUy meaningful curve is the 
blue solid one, which shows a completely different beha- 
viour. 

To the authors' knowledge, the only physical system 
that both couples directly to fermion fields and is easily 
measurable, the electromagnetic field, does so through 
an interaction quadratic on the ficld. The detector mod- 
els studied in [23, considered this kind of quadratic 
couplings. With some additional structure, these may 
well constitute physical models for a fermionic detector. 
However, these simplified detectors would not be able 
to acquire ali the entanglement from the fermionic field 
modes, save for a few very specific states. 

Of course, ali the considerations madc abovc apply 



to this fermionic detector as well as to any other, even 
though the detector is not capable of measuring field en- 
tanglement completely. 

V. CONCLUSIONS 

In this paper we tackle an ambiguity present in ali 
previous works on fermionic entanglement in non-inertial 
frames. The amount of entanglement is not invariant 
under nonlocal changes of basis. This is a general feature 
of quantum systems and hardly new. 

When working with fermionic systems as support for 
qubits, we must be carcful, since naively defining the 
basis of the multipartite Hilbert space without paying 
attention to the sign convention may lead to unphysical 
results. We bave shown that an entanglement measure 
can behave completely differently if we consider the dif- 
ferent bases naturally suggested by the fermionic oper- 
ator orderings. 

This situation has a special relevance in the context of 
fermionic fields in non-inertial frames, where part of the 
system has to bc traccd out. It could then be naturally 
expccted that once those modes bave been traced out, 
they would not affect entanglement of the remaining sys- 
tem. Nevertheless, we have seen that this intuition is not 
true. 

Of course, physics cannot depend on a mathematical 
convention. We have shown how when we introduce de- 
tectors, which are a fundamental requirement to observe 
the field, the ambiguity disappears. Field entanglement 
is only physical if it can be acknowledged by detectors 
coupled to the field. There is only one set of bases (ali of 
which share the same entanglement properties) for which 
field entanglement can be imprinted on detectors. We 
have shown the form of a simple detector model (consist- 
ing of an array of Unruh-Dewitt-like detectors) for which 
ali the entanglement can be transmittcd from the field to 
the detector. 

Incidentally, entanglement in this scenario happens to 
behave in a different way from that for the bases con- 
sidered in previous works in the literature on fermionic 
entanglement in non-inertial frames 0, H [llÌ . 
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